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ABSTRACT 

We use Luminous Red Galaxies from the Sloan Digital Sky Survey II to test the cosmologi- 
cal structure growth in two alternatives to the standard ACDM+GR cosmological model. We 
compare observed three-dimensional clustering in SDSS DR7 with theoretical predictions for 
the standard vanilla ACDM+GR model, Unified Dark Matter cosmologies and the normal 
branch DGR In computing the expected correlations in UDM cosmologies, we derive a pa- 
rameterized formula for the growth factor in these models. For our analysis we apply the 
methodology tested in Raccanelli et al. ( 20 1 Q) > and use the measurements of Samushia et al. 



( |201 that account for survey geometry, non-linear and wide-angle effects and the distribu- 
tion of pair orientation. We show that the estimate of the growth rate is potentially degenerate 
with wide-angle effects, meaning that extremely accurate measurements of the growth rate on 
large scales will need to take such effects into account. We use measurements of the zeroth 
and second order moments of the correlation function from SDSS DR7 data and the Large 
Suite of Dark Matter Simulations (LasDamas: McBri de et al.|20il) , and perform a likelihood 
analysis to constrain the parameters of the models. Using information on the clustering up to 
r max =120 Mpc/h, and after marginalizing over the bias, we find, for UDM models, a speed 
of sound ^ 6.1e-4, and, for the nDGP model, a cross-over scale r c ^ 340 Mpc, at 95% 
confidence level. 

Key words: large scale structure of the Universe — cosmological parameters — cosmology: 
observations — methods : analytical 



1 INTRODUCTION 

The strangest feature of our current cosmological model is the 
observation that the expansion rate of the Universe is accelerat- 
ing ( |Perlmutter et al.|1999"{|Riess et al.|1998) . Understanding the 
cause of cosmic acceleration is one of the great challenges of 
physics. It has been speculated that the cause of this acceleration 
is a cosmological constant, or perhaps some novel form of mat- 
ter; our ignorance is summarized by the simple name for the cause 
of the observed phenomenon: dark energy. Alternatively, it could 
be explained by the breakdown of Einstein's General Relativity 
(GR) theory of gravitation on cosmological scales (see |Durrer &| 
Maartens 2008 for a review on different dark energy and modi- 
fied gravity models). Observations of the large-scale structure of 
the Universe have played an important role in developing our stan- 
dard cosmological model and will play an essential role in our in- 
vestigations of the origin of cosmic acceleration. 

We will illustrate how it is possible to test GR and alternative 
models of gravity using Luminous Red Galaxies (LRG) from the 



Sloan Digital Sky Survey (SDSS) DR7 data. To estimate the sta- 
tistical errors on our measurements we use galaxy catalogues from 
the Large S uite of Dark Matter Simulations (LasDamas: |McBride] 
et al.|20lT[ lp] that are designed to model the clustering of Sloan 
Digital Sky Survey (SDSS) galaxies. 

The presence of a dark energy component in the energy- 
density of the Universe (or the fact that our theory of gravity needs 
to be modified on large scales), modifies the gravitational growth 
of large-scale structures. The large-scale structure we see traced by 
the distribution of galaxies arises through gravitational instability, 
which amplifies primordial fluctuations that originated in the very 
early Universe; the rate at which structure grows from small per- 
turbations offers a key discriminant between cosmological models, 
as different models predict measurable differences in the growth 
rate of large-scale structure with cosmic time (e.g. |Jain & Zhang 
[2007l[Song"& Koyama|2009||Song & Percival|2009| ). For instance? 



1 http://lss.phy.vanderbilt.edu/lasdamas/ 
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dark energy models in which general relativity is unmodified pre- 
dict different large-scale structure formation compared to Modified 
Gravity models with the same background expansion (e.g. Dvali, 
Gabadadze & Porrati|2000l |Carroll et al.|2004| [Brans et. al|2005| 
Nesseris & Perivolaropoulos 2008 , Ya mamoto et al.|2008|[20Tol >. 



Observations of Redshift Space Distortions (RSD) in spectro- 
scopic galaxy surveys are a promising way to study the pattern and 
the evolution of the Large Scale Structure of the Universe (Kaiser 
1987 , Hamilton 1998 ), as they provide constraints on the amplitude 
of peculiar velocities induced by structure growth, thereby allow- 
ing tests of the theory of gravity governing the growth of those 
perturbations. RSD have been measured using techniques based on 



Hawkins et al.|2003 


Percival et al.|2004||Pope et al.|2004||Zehav 


et al.||2005[ Tegmar 


k et al.||2006 Okumura et al.|2008| Cabr6 & 


Gaztanaga 2009 Guzzo et al.|2008; Blake et. al||201 1J>; the most 



[20121 [Sanchez et al.|2012) . 

A key element of RSD is that the motion of galaxies is inde- 
pendent of their properties and of the bias, that relates the bary- 
onic matter to the total mass; therefore, measurements of peculiar 
velocity directly probe the matter distribution. They also are com- 
plementary to other probes, since they depend on temporal metric 
perturbations, while e.g. weak lensing depends on the sum of the 
temporal and spatial metric perturbations and the Integrated Sachs- 
Wolfe effect depends on the sum of their derivatives. 

The standard analysis of RSD makes use of the so-called 
Kaiser formalism, that relies on some assumptions, including con- 
sidering only the linear regime and the distant observer approxi- 
mation, and restrains the range of usable scales to 30-60 Mpc/h. 
There have been several attempts to model smaller scales RSD, 
exploring the quasi-linear regime (e.g. Scoccimarro|2004| jTaruya, 



INishimichi & Saito|2010||Reid & White|2011||Kwan et al.|2011) ; 
recently Bertacca et ah] (20 12) developed a formalism to compute 
the correlation function including GR corrections, that arise when 
probing scales comparable to the Hubble scales. 

In this paper we show how precise measurements of the 
clustering of galaxies can be used to test cosmological models; 
we make use of the wide-angle methodology as tested in |Rac-| 
|canelli et al.| ( |2010) , that drops the distant observer approxima- 
tion, combined with prescriptions and measurements of SDSS-II 
data from Samu shia et al.|p011) , to constrain two interesting al- 
ternatives to the standard cosmology: a particular class of Unified 
Dark Matter (UDM) models (Bertacca et al.|2008| |20T0) and the 
normal-branch Dvali-Gabadadze-Porrati (DGP) (Schmidt 2009). In 
the process, we also derive a parameterized formula for the growth 
factor in UDM models, and we show that wide-angle corrections 
are degenerate with variations of the rate of the growth of struc- 
tures, demonstrating the need to include them if one wants to mea- 
sure the growth rate at percent level. 

The paper is organised as follows: in Section|2]we briefly re- 
view the theory of RSD; in Section [3] we present the catalog used; 
the methodology used to perform measurements is reviewed in Ap- 
pendix|A] in Section|4]we introduce the parameterisation of struc- 
ture growth we apply in our tests and discuss degeneracies that arise 
from a more comprehensive description of the data; after discussing 
the theoretical growth of structures for the non-standard cosmolog- 
ical models, we present the measurements in SectionB] and set lim- 
its on the parameters in Section|6] finally, in Section|7] we conclude 
and discuss our results. 



2 REDSHIFT-SPACE DISTORTIONS 

RSD arise because we infer galaxy distances from their redshifts 
using the Hubble law: the radial component of the peculiar velocity 
of individual galaxies will contribute to each redshift and will be 
misinterpreted as being cosmological in origin, thus altering our 
estimate of the distances to them. The relation between the redshift- 
space position s and real-space position r is: 



s(r) = r + v r {r)v , 



(1) 



where v r is the velocity in the radial direction. 

The measured clustering of galaxies will therefore be 
anisotropic and the additional radial signal can be used to deter- 
mine the characteristic amplitude of the pair-wise distribution of 
the peculiar velocities at a given scale, which in turn depends on 
the growth rate. 

Measurements are normally obtained over a small range of 
scales, because of simplified modeling. In this work, we use the ex- 
tended analysis tested in Raccan elli et al.|(2010) and |Samushia et| 
[aL](201 1) , which includes a more realistic description of the geom- 
etry of the system, dropping the plane-parallel approximation; this 
allows us to fit the observed galaxy correlation function on a larger 
range of scales, and therefore to be more sensitive to the cosmolog- 
ical parameter variations. 

By imposing the conservation of the number of galaxies we 
can derive the lacobian for the real- to redshift-space transforma- 
tion (at the linear order): 



<T(s) = *"(r) - 



dv a(r)v 
dr r 



(2) 



where S s ' r are the observed redshift- and real-space galaxy over- 
density at positions s and r, and: 

dlnr 2 N r (r) 
dlnr ' 

and N r (r) is the expected galaxy distribution in real space. The 
simplest statistic that can be constructed from the overdensity field 
is the correlation function { (ri2), defined as: 



x(p) = 



(3) 



«na) = <*(pi)*(pa)>. 



(4) 



In linear theory, all of the information is enclosed in the first three 
even coefficients of the Legendre polynomial expansion of the 
function £ (Hamilton|1992) : 

£(r, /0 = &(r)Lo(/*) + £ 2 (r)L 2 (p) + Ur)Mn) , (5) 

where Lt are the Legendre polynomials and y, is the cosine of the 
angle with the line of sight (jj, = cos(4>) in Figure All. 



In this work we will use measurements of the momenta of 
the correlation function from the SDSS DR7 catalogue, using a 
methodology that takes care of several corrections, as we will de- 
scribe in the next Sections. 



3 THE SLOAN DIGITAL SKY SURVEY 

We use data from the SDSS-II data release 7 (DR7), which obtained 
wide-field CCD photometry (Gunn et al.| 1 9 98 ) in five passbands 
(u, g,r,i, z; e.g. Fukugita et al. 1996 1, amassing nearly 10,000 
square degrees of imaging data for which the object detection is 
reliable down to r ~ 22 ( Abazajian et al. 2009). From these photo- 
metric data, Luminous Red Galaxies (LRG) were targeted (Eisen-| 
|stein et al.|2001) and spectroscopically observed, yielding a sample 
of 106,341 LRGs in the redshift bin 0.16 < z < 0.44. 
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An estimate of the statistical errors associated with the mea- 
surements is achieved through LasDamas mock catalogues, which 
model the clustering of Sloan Digital Sky Survey (SDSS) galaxies 
in the redshift span 0.16 < z < 0.44. The simulations are pro- 
duced by placing artificial galaxies inside dark matter halos using 
an HOD with parameters measured from the SDSS galaxy sample. 
We use the 80 "Oriana" catalogs that have exactly the same angular 
mask as the SDSS survey and subsample them to match the redshift 
distribution of the Luminous Red Galaxies in our SDSS DR7 data 
set. Further details regarding the angular and redshift distribution 
of galaxies in random catalogues can be found in Samushia et al. 
pOTT) . 

3.1 Methodology 

To perform our analysis of cosmological models we use the 
methodology presented in |Raccanelli et aL 1 2010} a nd measure- 
ments of SDSS-II data from Samushia et al. l 2011| l. In the Ap- 
pendix we briefly revisit the main aspects of the approach we fol- 
lowed; this drops the distant observer approximation and includes 
a careful treatment of survey geometry, non-linear effects and the 
distribution of pair orientation. This allows us to consider a wide 
range of scales (we use measurements from 30 to 120 Mpc/h). 



4 PARAMETERIZING THE GROWTH OF STRUCTURE 

Measuring the matter velocity field at the locations of the galaxies 
gives an unbiased measurement of f(Ts m , provided that the distri- 
bution of galaxies randomly samples matter velocities, where: 



f 



dlnD 
d In a 



(6) 



is the logarithmic derivative of the linear growth rate, D(a) oc S m , 
with respect to the scale factor a (5 m being the fractional matter 
density perturbation) and ag m quantifies the amplitude of fluctua- 
tions in the matter density field. 

|Linder| ( [2005| > proposed a gravitational growth rate formalism, 
which parameterises the growth factor as: 



D(a) — a exp 



[n^(a') 



da' 



which leads to the following expression for /: 



with: 



/=[n m (o)] 



E^exp^MaO + l]^ 



(7) 



(8) 



(9) 



where the summation index goes over all the components of 
the Universe (i.e. dark matter, dark energy, curvature, radiation). 
Within this formalism, 7 is a parameter that is different for different 
cosmological models: in the standard ACDM+GR model it is con- 
stant, 7 as 0.55, while it is ~ 0.68 for the self-accelerating DGP 
model (see e.g. finder 20051. In some other cases, it is a function 
of the cosmological parameters or redshift, as we will discuss. It 
should be noted, however, that the parameterization given by Equa- 
tion|7]does not necessarily describe the growth rate in non-standard 
cosmologies. 

Given that, as we will see, depend on /, measuring RSD 
allows us to determine 7, and hence to test different cosmological 
models. This also provides a good discriminant between modified 




Figure 1. Mock 6 distributions used to illustrate effects of neglecting wide- 
angle corrections when measuring 7 . 



gravity and dark energy models, as argued by Linder (2005 2007) ); 
|Guzzo et aLlpOQg) . 



4.1 Degeneracy #-7 

In Appendix |A| we review the methodology used to obtain mea- 
surements of multipoles of the correlation function; our method- 
ology involves corrections with respect to standard analyses, and 
in particular it drops the distant observer approximation (6 = 
in Figure |A~T) , The error in the estimate of the correlation function 
induced by assuming 8 = can lead to a wrong estimate of the cos- 
mological parameters measured with Redshift-Space Distortions. 

We investigate the systematics introduced by assuming 9 = 
for four different mock distributions of 6 (Dl, D2, D3, _D4). In 
Figure [T] we plot the galaxy distribution as a function of the sepa- 
ration angle for each of them. The ratio between the monopole of 
the correlation function, £n> derived under the Kaiser approxima- 
tion and those obtained when including the wide-angle correction 
is shown in Figure[2] as a function of separation scale for the stan- 
dard cosmological model. The same quantity for the quadupole of 
the correlation function, £2, is plotted in Figure|5] As expected, the 
largest errors are induced by the distributions D2 and D4, which 
have galaxy pair distributions peaked at relatively large 8; however, 
we can note that even the D 1 distribution introduces a correction of 
the order of 5%. For comparison, we also add the effect that a vari- 
ation of the 7 parameter would have in the Kaiser approximation. 
Even though the shape at large angular separation is different, the 
size of the variation is comparable and the wide-angle correction 
must be included for a proper, more exact, analysis. 

The distributions used here are illustrative; a more careful 
analysis of the influence of wide-angle and other large-scale cor- 
rections for future Euclid-like and SKA-like surveys is left to a 
follow-up paper. 



5 MEASUREMENTS 

The momenta of the correlation function are sensitive to the 7 pa- 
rameter through the function /. In this work we concentrate on 
two alternatives to the standard ACDM+GR scenario: i) the Unified 
Dark Matter (UDM) cosmology ( |Bertacca et al.|2008||2010| >, and 
ii) the normal branch DGP model including Dark Energy (nDGP) 
of [Schmidt (2009). These two models deviate from the standard 
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Figure 2. Ratio of monopole of the correlation function in the ACDM+GR 
model in the Kaiser approximation to the monopole of the correlation func- 
tion including wide-angle corrections for the mock (^-distributions of Fig- 
ure^ compared with the effect of varying 7 in the Kaiser approximation. 
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Figure 4. Wide-angle corrections (i.e. £o( r i 9 )/fo( r i 6 = 0)) to the 
monopole of the correlation function for the models considered and dif- 
ferent values of their parameter. 
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Figure 3. Ratio of quadrupole of the correlation function in the ACDM+GR 
model in the Kaiser approximation to the quadrupole of the correlation 
function including wide-angle corrections for the mock ©-distributions of 
Figure [T] compared with the effect of varying 7 in the Kaiser approxima- 
tion. 



cosmology in different ways: the UDM model assumes a single 
dark fluid with a clustering part, and leaves GR unmodified, while 
in nDGP, gravity crosses over from 4D to 5D above a cross-over 
scale r c (see Sections [5.1| and |5^2| for details). For the UDM model, 
we use the values of 7 of Equation [T3] for a set of values of the 
speed of sound , while for nDGP we compute the growth solv- 
ing explicitly Equation|17| 

As one can see in Equation |A4| the wide-angle and mode- 
coupling corrections are described by a set of terms that depend 
on {r, 6, (j>, 7}, and this means that wide-angle corrections are also 
different for different models of gravity. In Figure [4] [5] we show 
the wide-angle corrections to the monopole and quadrupole of 
the correlation function for the different models, compared to the 
ACDM+GR case. 



5.1 UDM 

Assuming a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) 
background metric with scale factor a(t), Bertacca et al. j2008 
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Figure 5. Wide-angle corrections (i.e. £2 (r, 6)/(, 2 (r, 9 = 0)) to the 
quadrupole of the correlation function for the models considered and dif- 
ferent values of their parameter. 



|2010[ >, introduced a class of Unified Dark Matter (UDM) scalar 
field models which, by allowing a pressure equal to — c 2 pA on cos- 
mological scales, reproduces the same background expansion as 
the ACDM one. In such a way, a single scalar field can be respon- 
sible for both the late time accelerated expansion for the Universe 
and of the growth of structure. When the energy density of radia- 
tion becomes negligible, the background evolution of the Universe 
is completely described by: 



H 2 {z)=Hl [« A0 + OnO (1 + Zf] 



(10) 



where H is the Hubble parameter, Hq = H(z = 0) and z is the 
redshift. Qao an d QmO = 1 — ^ao can be interpreted as the "cos- 
mological constant" and "dark matter" density parameters, respec- 
tively. 

The density contrast of the clustering component is <5dm = 
Sp/poM, where pdm = P — Pa is the only component of the 
scalar field density that clusters. In these models one of the most 
relevant parameters is the sound speed of the perturbations, which 
defines a typical sound-horizon (Jeans length) scale above which 
growth of structure is possible (Bertacca & Bartolo 2007 1. Follow- 
ing Bertacca et al. < 20 1 IJ, for scales smaller than the cosmological 
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Figure 6. Deviation of the growth factor for the UDM model, as a function 
of the speed of sound Coo , for different values of k. 



Figure 7. 7tjdm as a function of k, for different values of c a 
0.15), compared with the ACDM+GR case. 



(at 2 



horizon and z < z ICC , we have that: 

5dm [k; T](z)] = Tudm [fc; 77(2)] S m [k; r)(z)] 



(11) 



where 8 m is the matter density perturbation in the standard ACDM 
model, rj is the conformal time and Tudm(&; rj) is the transfer func- 
tion for the UDM model: 



Tudm(M) =jo[A{ri)k] 



1(a) = 



Cs(r]') dr)', 



+ (1 — c^)fi m oa 



(12) 
(13) 

(14) 



We then define the parameter Coo as the value of the sound speed 
when a — > 00. 

In this model, the growth factor is computed correcting the 
ACDM one with the transfer function of Equation[T2] its deviation 
from the standard model case, as a function of Coo for different 
values of k, is plotted in FigurelS] 

We refer to |Bertacca eTaflp008| [20T0] |20Tfl > for details of 
the UDM models considered (see also |Camera et al.|2009| |2010| 
|Piattella & Bertacca|20lT1|Camera et al.|2012) . 

Here we report a new result, which is particularly relevant 
for the present analysis: the parameterization of the growth rate in 
UDM models, which turns out to be: 



In 



7udm(o., k, c 



dlnT UDM (q,fc,c 
d In a 



+ J*ACDM 



lnfi m (a) 



(15) 



and in Figure [7] is 



as one can see, 7udm depends on (a, k, c a 
shown its value as a function of them. 

Equation[T5]gives the value of 7 for UDM cosmologies; it de- 
pends on redshift, scale and speed of sound. Figure [7] shows that 
7udm 7acdm- Actually this is easy to understand: one of the 
features of the UDM models is that, under the Jeans length, the den- 
sity contrast decreases while oscillating in time (see, e.g., |Bertacca| 
|& Bartolo|2007l|Bertacca et al.|2011J . This means that, compared 
to the ACDM model, there is a further suppression in the growth of 
structures. Notice also that there are some values of Coo that give 
non-physical results when one uses the parameterization of Equa- 
tion [T5] however, even using Equation [6] the growth rate for these 
values becomes highly oscillatory, so we will consider them ruled 
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Figure 8. £0 measured from SDSS DR7 and theoretical predictions for 
ACDM and UDM with different values of the speed of sound. 



out. For this reason, given that we will limit our analysis to lin- 
ear scales, and accounting also for existing limits on UDM models 
(e.g. |Bertacca et al.|201 1| > we will consider values of Coo < 0.002. 

In Figures |8|9| we show the comparison of measurements of 
the zeroth and second order momenta of the correlation function 
measured from SDSS DR7 with the theoretical predictions for 
UDM models, for several values of the speed of sound parameter 

Coo • 

5.2 DGP 



In the Dvali-Gabadadze-Porrati (DGP, |Dvali, Ga badadze & Por-| 
|rati||200 0) model, all matter and radiation are confined to a four- 
dimensional brane in five-dimensional Minkowski space. Gravity, 
while restricted to the brane on small scales, propagates into the 
extra dimension above the cross-over scale r c . This scenario ad- 
mits an FRW cosmology on the brane (Deffayet 2001 1, where the 
Friedmann equation is modified to: 



u2 H 87TG r _ 



+ PdeJ 



(16) 



The sign on the left-hand-side depends on the choice of embedding 
of the brane. The negative sign (accelerating branch) leads to an 
accelerated expansion of the Universe at late times without a cos- 
mological constant or Dark Energy (Deffayet 200TJ, i.e. pDE = 
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Figure 9. £2 measured from SDSS DR7 and theoretical predictions for 
ACDM and UDM with different values of the speed of sound. 



Figure 10. Deviation of the growth factor for the nDGP model, as a function 
of the cross-over scale r c (see text for details). 



and r c ~ Hg 1 , while the positive sign (normal branch) does not 
yield acceleration by itself. The simplest self-accelerating model 
is in conflict with observations of the CMB and Supernovae (e.g., 
Lomb riser et al.|2 009 1, and also has theoretical issues ( |Luty et aT] 
2003[|Nicolis & Rattazzip004{ [Gregory et al.|2007| >. 

Here we consider a normal-branch DGP model including a 
Dark Energy component pde to yield an accelerated expansion. 
Specifically, we use the model of Schmidt ( 2009 1 where the equa- 
tion of state of the Dark Energy is tailored to yield an expansion 
history identical to ACDM for all r c . This model is not ruled out 
by expansion history probes or theoretical issues. Furthermore, it 
can serve as a toy model for more recent scenarios ( deRham et al.| 
[2008l|Afshordi et al.|2008) . 

On scales much smaller than the horizon and smaller than r c , 
but yet large enough so that linear perturbation theory applies, the 
growth of perturbations during matter domination in these models 
is governed by (Koyama & Maartens 2006): 

S + 2H5 = ^n m (a)a 2 H 2 (l + -^—\s, (17) 

where dots denote derivatives with respect to time, and 

/3 = l±2Hr c (l + ^p). (18) 

Here the positive (negative) sign holds for the normal (self- 
accelerating) branch, respectively. Note that in the latter case, /3 < 
and hence G c g < G, i.e. the growth of structure is slowed down 
with respect to GR. The opposite is the case for the normal branch 
(/? > 0). In the case of the self-accelerating (sDGP) model with 
Pde — 0, Linder ( 2005 ) showed that the growth rate can be well 
described by the parameterization / s dgp(cs) = f2 m (a) 7sDGP , with 
7sDGP = 0.68. For the normal branch model with Dark Energy 
(nDGP), this parametrization does not provide a good fit, and we 
use a numerical integration of Equation [T7] to derive the growth 
rate. The change in the growth with respect to ACDM is shown in 
Figure 10 as a function of r c . 

In Figures [TT][T2] we show the comparison of the zeroth and 
second order momenta of the correlation function measured from 
SDSS DR7 with the theoretical prediction for the nDGP model, for 
some values of the cross-over scale r c . 
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Figure 11. £0 measured from SDSS DR7 and theoretical predictions for 
ACDM and nDGP, for some values of the cross-over scale r c . 

6 RESULTS 

We compute the likelihood of the and 1 /r c parameter for UDM 
and nDGP model respectively, conditioned to the other parameters 
fixed to the WMAP 7-year best-fit ACDM model ones. We as- 
sume a Gaussian likelihood with covariance matrix described by 
Equation |A10| This choice has the advantage of having the ACDM 
model as an asymptotic limit, when the extra parameter tends to 
zero. 

We evaluated a joint C for the zeroth and second momenta 
of the correlation function, focusing on two cases: (i) r max = 80 
Mpc/h, where the data start to deviate from the mock for both £0 
and £2 (see Figure [A4]|A5[ l, and (ii) r max = 120 Mpc/h, in order to 
investigate larger scales and where the data are still in reasonable 
agreement with the mocks. For larger scales, Ross et al. ( 201 1 ) sug- 
gested (when considering a different sample) that the excess power 
can be due to systematic errors, so, to be conservative, we will not 
fit these scales. 

The bias is poorly understood: in order to take into account 
its theoretical uncertainty, we leave it as a free parameter and 
marginalize over it. Since additional information on the bias can 
be independently inferred from other probes, as for example lens- 
ing measurements, we also derive the constraints on cosmological 
parameters which can be obtained when fixing the bias at its best-fit 
value assuming a ACDM cosmology. 



© 0000 RAS, MNRAS 000, 000-000 



Testing Gravity using RSD 7 




40 60 80 100 120 1 2 3 4 5><10" 

r (Mpc\h) l/r c (Mpc) 



Figure 12. £2 measured from SDSS DR7 and theoretical predictions for 
ACDM and nDGP, for some values of the cross-over scale r c . 



Figure 14. Likelihood for the nDGP model, as function of the inverse of the 
cross-over scale l/r c , when assuming knowledge of the bias (solid line), 
and after marginalizing over it (dashed line). 



— fixed bias 




Figure 13. Likelihood for the UDM model, as function of the speed of 
sound Coo , when assuming knowledge of the bias (solid line), and after 
marginalizing over it (dashed line). 



In Figure s [T3| and [T4| we show the likelihood, as a function of 
Coo (UDM) and l/r c (nDGP), in the r max = 120 Mpc/h case, as- 
suming the knowledge of the bias (solid lines), and after marginal- 
izing over it (dashed lines). 

In Table[T]we report the best fit and the constraints on Coo and 
r c at different confidence levels, for the various cases considered. 

The bias plays a very important role: this was expected, since 
the shape of the correlation function is a weak function of the addi- 
tional parameter in both classes of models, whereas the amplitude 
variation is large (in UDM the shape variation is more relevant, and 
increasingly important for larger values of Coo )■ 

In the UDM case, assuming the knowledge of the bias re- 
sults in a best fit value of Coo = 0, which corresponds to the 
ACDM model, while, after marginalizing over the bias, the best fit 
is Coo = 5 x 1CP 4 . From the Table [T] emerges that the constraints 
on the parameters are not sensitive to the maximum value of the 
distance between galaxy pairs. This is not surprising, as the devi- 
ation of the model from ACDM is larger for smaller scales (this 
can also be seen from Figures[8][9|. It is interesting to note that, af- 
ter bias marginalisation, the best fit does not correspond to ACDM 
, meaning that the shape of the measured correlation functions is 
not fitted by ACDM ; it would be interesting to perform a future 
analysis at different redshift or with other data sets (e.g. BOSS), to 



Model 


l-o- 


2-cr 


3-cr 


Best Fit 


UDM[80], 


Coo «S 2.1e-4 


Coo s£ 3.8e-4 


Coo «S 4.7e-4 


Coo = 


fixed bias 








(ACDM ) 


UDM[80] 


Coo ^ 5.8e-4 


Coo «S 6.3e-4 


Coo «S 6.9e-4 


Coo = 5.0e-4 


marginalized 










UDM[120], 


Coo «S l-9e-4 


Coo «S 3.5e-4 


Coo «S 4.3e-4 


Coo = 


fixed bias 








(ACDM )4 


UDM[120] 


Coo 5.8e-4 


Coo ^ 6.1e-4 


Coo ^ 6.9e-4 


Coo = 5.0e-4 


marginalized 










nDGP[80], 


r c > 12580 


r c > 7040 


r c 5> 5500 


r c = 25000 


fixed bias 


(Mpc) 


(Mpc) 


(Mpc) 


(Mpc) 


nDGP[80] 


r c > 660 


r c ^ 290 


r c ^ 255 


r c = 00 


marginalized 


(Mpc) 


(Mpc) 


(Mpc) 


(ACDM ) 


nDGP[120], 


r c > 9803 


r c > 6480 


r c ^ 5050 


r c = 20000 


fixed bias 


(Mpc) 


(Mpc) 


(Mpc) 


(Mpc) 


nDGP[120] 


r c ^ 1237 


r c > 340 


r c ^ 270 


r c = 00 


marginalized 


(Mpc) 


(Mpc) 


(Mpc) 


(ACDM ) 



Table 1. Constraints and best fits for models considered, when knowledge 
of the bias is assumed and after marginalizing over it as a free parameter, 
for the r ma x=80 and 120 Mpc/h cases. 



see how the peak of the likelihood will be modified. Our final con- 
straint at 95% confidence level, after marginalizing over the bias, 
is Coo ^ 6.1 X 10~ 4 , almost two orders of magnitude better than 
previous constraints (Bert acca et al.|20TT) . 

When considering the nDGP scenario, the effect of marginal- 
izing over the bias is even more dramatic: the constraints on r c 
worsen by a factor of ~20 (see Figure [14] and Table [TJ. This oc- 
curs because the deviation from ACDM+GR of this model is scale 
independent (below the cross-over scale), hence degenerate with 
the bias. After bias marginalisation, the best fit corresponds to the 
ACDM+GR model. However, for this model there are no published 
constraints in the literature, so our analysis presents a first result on 
that. Note that a cross-over scale of r c ~ 300 Mpc implies strong 
modifications to gravity on larger scales, and the integrated Sachs- 
Wolfe (ISW) effect in the cosmic microwave background is likely 
able to constrain such values as well. 
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7 CONCLUSIONS 

In this work, we showed that we can use Redshift-Space Distor- 
tions to test cosmological models, by measuring the monopole and 
quadrupole of the correlation function of galaxies, {^0,^2}- The 
methodology and robust measurement of the correlation function 
presented in Raccanelli et al. (20101; Samushia et al. (2011), which 
includes a careful treatment of corrections due to the geometry of 
the system, allows us to use those multipoles in a wider range of 
scales w.r.t. most standard analyses, and so have better constraints 
on the models tested. As shown in Samus hia et al.| l |2011) , most 
of the approximations assumed in the Kaiser analysis need to be 
dropped to have precise measurements of the correlation function, 
and so of the growth rate, whose deviation from the ACDM+GR 
expected value would imply the need of a new cosmological model. 

We explored the Unified Dark Matter (Bertacca et al. 12008) 
|2010| > and the normal branch DGP ( Schmidt 2009|> models, that 
present deviations from the ACDM+GR scenario. Both classes of 
models are parameterised by one additional number: the speed of 
sound Coo and the cross-over scale r c , for UDM and nDGP respec- 
tively. The value of these parameters affects both the growth rate 
parameter 7 and the wide-angle corrections. Moreover, UDM mod- 
els are characterised by a growth rate parameter, 7, that depends 
on {k, z, Coo}. After deriving its analytic expression, we used it to 
compute the predicted momenta of the correlation function. 

We then compared observations of LRGs from SDSS-II DR7 
with theoretical predictions from the two cosmologies. This anal- 
ysis allowed us to tighten the constraints on the speed of sound of 
UDM models to Coo ^ 6.1 x 10~ 4 , and to put, for the first time, 
a lower bound on the cross-over scale for the nDGP model of 340 
Mpc, both at 95% confidence level. It is worth noting that the re- 
sults would largely benefit from a better knowledge of the bias that 
could be obtained combining information from other cosmological 
probes, as e.g. gravitational lensing. 

We showed the potential of this methodology for constraining 
alternative cosmological models, in particular when future surveys 
will provide access to a wider range of scales, hence allowing a 
tomographic analysis, which will be essential to break degeneracy 
between cosmological parameters. 
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APPENDIX A: METHODOLOGY 
Al The wide-angle corrections 

Most previous RSD analyses have used the simple plane-parallel 
approximation given by the Kaiser formula. In this case, a Fourier 



mode 5 s (k) in redshift space is simply equal to the unredshifted 
mode <5(k) amplified by a factor 1 + f3fJ%. : 



6 s (k) = (l + p£)5(k), 



(Al) 



where (3 = f/b, with 6 being the bias. This correction arises from 
the Jacobian of Equation^ when the (1 + term is neglected. 
The wide-angle linear redshift-space correlation function and 



power spectrum have been analytically derived by |Zaroubi & 
|Hoffman|([T993l>; |Szalay etal] p998) ; |Szapudi| pOOl) - JMatsubara 



( 2004lTfPapai & Szapudi] (|2008| >, and tested against both simula 
tions (Raccanelli et al. 2010) and real data ( |Samushia et al.|201 1) . 

Papai & Szapudi (2008) have argued that, for wide angles, 
the v/r term in Equation [2] is of the same order as the d r v term. 
As a consequence, a careful and precise analysis of the correlation 
function, requires the full (linear) Jacobian, dropping the distant- 
observer approximation. In this case, we can express the linear 
overdensity as: 



5 s (s) 



e 3 3 



(2tt 



1 + m • k, 



ia(r)f 



*3 ' kj 

rk 



6(k); 
(A2) 



the redshift-space correlation function reads: 



r = <<5 s ( Sl )5 s *(s 2 )) 
/ , 2/ 



(Tk 
(2tt) 3 



P(k)e 



ifc(ri — r 2 ) 



o r 2 k 



(A3) 



The third terms in the brackets describe the wide-angle effects, 
while the fourth ones are responsible for the so called mode- 
coupling. The r\ and r 2 terms in the denominator depend on the 
angular separation of the galaxies, and a is proportional to the loga- 
rithmic derivatives of the galaxy distribution function (Equation^. 
We refer to Figure [AT] for the geometry of the problem, where we 
define 20 to be the angular separation of the two galaxies consid- 
ered, (j>i as the angle between the vector to the first galaxy in a pair 
ri and r, r to be the vector connecting galaxies in a pair, and <f> 2 to 
be the angle between vector to the second galaxy in a pair r2 and 
r. 

Tripolar spherical harmonics are the most natural basis for 



the expansion of a function that depends on three directions ( |Var- 
shalovichetal.11988), so, as suggested by Szapudi ( 2004 1 and Papai 



& Szapudi ( 2008 ), we expand Equation |A3| using a subset of them, 



so that the redshift-space correlation function can then be written 

as: 



e 2 e(ri,r 2 ,r) 



(A4) 



where B eil2l (r, cj> u <p 2 ) arc a scries of coefficients that depend on 
/, gi(4>i) an d £H r ) (see |Raccanelli et al.|2010 for details on the 
definition of these functions). These coefficients can be divided into 
two different subsets: one that depends on the third term inside the 
brackets in Equation 



A3 



given by B 1 2 (r), with l\l 2 l combi- 
nations of 0, 2, 4, with only the radial dependence accounting for 
the wide angle effects, and o ne th at depends on the fourth terms in- 
side the brackets of Equation A3 given by S^(r,0i,?i.2),with 
£i£ 2 £ combinations of 0, 1, 2, 3, with also an angular dependence, 



describing the mode coupling part (see Sz alay et al.| ( |1998"] >, |Sza-| 
|pudiH2004) and |Papai & SzapudiH2008) for a detailed derivation). 
The plane-parallel approximation emerges as a limit when fi = r 2 . 
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O 

Figure Al. The coordinate system adopted for the triangle formed by the 
observer O, and galaxies Gi and G2- 

This formalism was shown to accurately reproduce wide-angle ef- 
fects seen in numerical simulations (Raccanelli e t al.|2010| >. 

A2 SDSS DR7 

An extensive analysis of the SDSS DR7 data is provided in 
|Samushia et al.| ( [201 1[ >; we refer the reader to that paper for the 
full details on the analysis of SDSS DR7 data and for the relative 
importance of various corrections considered. In the following we 
use the measurements of momenta of the correlation function pre- 
sented there, to test models of gravity. 
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Figure A2. Measured distribution of angular separation 9 for LRG pairs in 
the SDSS DR7 catalog, at different scales. 
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Figure A3. Measured distribution of orientation angle fj, for LRG pairs in 
the SDSS DR7 catalog, at different scales. 



can be applied to both plane-parallel and wide-angle analyses. In 
Figure [A3] is shown the distribution of fi for observed LRGs in the 
SDSS DR7 catalog. 



A2.2 Non-linear Effects 

We model two non-linear effects: one due to the Baryonic Acoustic 
Oscillations (BAO) peak and the other due to the so-called Fingers 
of God. Since we are only interested in the signal on large-scales, 
where the linear theory is an accurate description, we assume that 
the non-linear effects are small, except for the fact that the matter 
power spectrum itself goes non-linear. We describe the non-linear 
contribution of the power spectrum by means of a two-component 
model, which splits P(k) into a "smooth" part that describes the 
overall shape and a "wiggled" part, which describes the BAO: 

fBAC>(fc, fi) = ffull(fc, M) - -Psmooth(fc, At) ■ (A5) 

The primary non-linear effect on the BAO component of the 
power spectrum is a damping on small scales, which can be well 
approximated by a Gaussian smoothing dBharadwaj 1996, Crocce 



A2.1 Distribution of 9 and fi 

For surveys that cover a significant fraction of the sky, the distri- 
bution of galaxies pairs has a complicated dependence on the vari- 
ables {r, fi, 9}, since not all sets of their combinations are equally 
likely or even geometrically possible. In particular, the distribution 
of fi does not correspond to that of an isotropic pair distribution, 
and this will strongly bias measurements of angular momenta of 
the correlation function. We can see in SectionfA|the effects on the 
correlation function of a non-zero fixed angular separation, and in 
Section pnj that these errors can bias our estimate of the 7 parame- 
ter. 

In Figure [A2] is shown the distribution of 9 of observed LRGs 
in the SDSS DR7 catalog. As expected, in this case the number of 
pairs with relatively large 9 is very small, and it increases with the 
linear separation of pairs. 

Other than the wide-angle and the mode coupling corrections, 
there is also another difference that has to be taken in account when 
doing a real-data wide-angle analysis, that derives from the fact 
that the distribution of galaxies in fi will be non-trivial, with some 
values of fi not permitted for non-zero 9. As a consequence, we will 
not be able to measure pure Legendre momenta of the correlation 
function, but instead we will need to use weighted integrals and 
biased momenta; corrections due to a non-uniform "/i-distribution" 



|& Scoccimarrol2 006 2008 1 |Eisenstein, Sep & W hite 2007j |Mat- 
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|subara|2008a|bl|Matarrese & Pietroni|2007l >: 



-PbaoO, (A = p b£o(k, m) ex P \ ~ k 



+ 



2 2 

(A6) 

where Ex = EnD and En = Eo(l + f)D; Eo is a constant 
phenomenologically describing the diffusion of the BAO peak due 
to non-linear evolution. From N-body simulations its numerical 
value is of order 10 Mpc/h and seems to depend linearly on erg but 
only weakly on k and other cosmological parameters. 

Within dark matter haloes the peculiar velocities of galaxies 
are highly non-linear. These velocities can induce RSD that are 
larger than the real-space distance between galaxies within the halo. 
This gives rise to the observed fingers of god (FOG) effect, that 
is a strong elongation of structures along the line of sight ( |Jack-| 
|son|1972) . The FOG effect sharply reduces the power spectrum on 
small scales compared to the predictions of the linear model, and 
is usually modeled by multiplying the linear power-spectrum by 
a function F, that depends on the average velocity dispersion of 
galaxies within the relevant haloes, cr v , k and fi. The most common 
one is a Gaussian damping (see e.g. |Cole et al.| 1995 Peac ock &| 
|Dodds|1996| >: 



F(a v , k, fx) = exp [— (ka v fx) 2 ] 



(AT) 



that is small on small scales and approaches unity for scales larger 
than l/a v . 

A3 The Estimator of the Correlation Function Momenta 

To estimate the momenta of the correlation function we use Landy- 
Szalay type estimators (Landy & Szafa yfl993} : 



E 3 , fc [DD{ n ,ii h 



(A8) 

2DR(n, fij , k ) + RR(n, m , 6 k )] 



where fx=cos(<p), while DD(n, fij,9k), DR(ri, fXj ,8k) and 
RR(ri, , 9k) are the number of galaxy-galaxy, galaxy-random 
and random-random pairs in bins centered on n, Lij and 8 k- 

RSD measurements are often extracted from the normalised 
quadrupole Q (Hamilton 1992), defined as: 

Ur) 



Q(r) 



6>(r) 



io{r')r' 2 dr' 



(A9) 



The normalised quadrupole Q was introduced because it is inde- 
pendent of the shape of the power spectrum, and so it depends only 
on the ft parameters, allowing to directly test gravity; however, this 
is true only in the Kaiser analysis, and so it is not true in our case. 
For this reason we test different cosmological models fitting the 
momenta of the correlation function, £o, £2- 

We compute error bars as the square root of the diagonal terms 
in the covariance matrix: 

C = ^ E - X M - > ( A10 > 

where X(r) is a vector of the measurements of £q , £2 at scale r and 
X is the mean value from all 80 mock catalogues. 

We compare measurements with predictions from different 
models, where we compute the redshift-space correlation function 
including effects from wide-angle and /i-distribution as well as sur- 
vey geometry and non-linearities. Estimates of Legendre momenta 



given by Equation |A8| correspond to: 

6(r)= f W I {r,9,^)i s {r,9,y)L l ^)d9d^ (All) 



where W r {r, 9, li) is a weight function that appears because of the 
geometrical constraints and the not uniform distribution of 11, as 
explained in Section [AT] and gives the relative number of pairs in 
a survey that form angles 11 and 9 for a given scale r; when the 8 
distribution tends towards a delta function centred at = 0, the 
wide-angle effects become negligible, while when the distribution 
in fx tends towards a uniform one, this effect becomes negligible. 
The correlation function of Equation |A4| can be written as: 

£ s {r,e,<p) =^c ab {f,r)L a \cos{6)]L h [cos{ V )l (A12) 



where {r, 9, ip} are as in Fig. Al and c a b(f, r) are coefficients that 



depend on the separation between galaxies and the / of Equation[6 



(see[Szapudi 2004{ |Papai & Szapudi|2008"l |Raccanelli et al.|2010 
for further details). 



A4 Las Damas Mocks 

In this paragraph we show our measurements from the LasDamas 
simulations and the ACDM+GR model prediction. As one can see 
in Figures |A4| |A5| the LasDamas set reproduces quite well the 
measured momenta of the correlation function, the main deviation 
being on large scales of the monopole of the correlation function, 
that represents also the main source of deviation from ACDM+GR 
at the scales we considered; however the other models don't help 
much in picking up that deviation. This is a well known discrep- 
ancy between theoretical predictions and observations, and it has 
been detected in spectr oscopic JKazin et al.|2010 Samush ia et al.| 
|201 1\ and photometric ^Thomas et al.|2010 1 data sets. This excess 
power at large scale can be induced by primordial non-Gaussianity 
(see e .g. |Matarrese & Verde|2008| |Dalal et al.|2008) |Slosar et al"] 
2008; Desjacques & Seljak 2010, Xia et al. 2010) or exotic physics 
(Thomas et al. 2010); however Samushia et al. 2011 found that 
the redshift dependence of the excess power is different from what 
would be caused by non-Gaussianity, and Ross et al. ( 20 1 1 1 suggest 
that this is likely to be due to masking effects from stellar sources. 
After correcting for systematics, Ross et alT] {2012| ) found consis- 
tency at better than 2a between the BOSS CMASS DR9 large-scale 
clustering data and the WMAP LCDM cosmological model. Fur- 
ther tests using mocks suggested that there was no evidence that 
additional potential systematic trends contributed individually at a 
level above that expected through noise. But, there is always the 
potential for more systematic problems to be present, or for the 
combination (which always tends to add power) not to have a more- 
significant contribution - i.e. form part of the 2a discrepancy. 

Figures |A6| |A7| show that our methodology, that includes 
corrections due to non-linearities, wide-angle and /i-distributions, 
can fit the ACDM+GR simulations very well, up to scales of 180 
Mpc/h; For the data analyses, however, we decided to use as a max- 
imum scale r = 120 Mpc/h. 



REFERENCES 

Abazajian, K. N., et al., ApJS, 182, 543, 2009 

Afshordi, N. and Geshnizjani, G. and Khoury, J. arXiv:astro- 

ph/08 12.2244, 2008 
Bertacca D., Bartolo N., JCAP, 071 1, 026, 2007 



© 0000 RAS, MNRAS 000, 000-000 



Testing Gravity using RSD 1 1 




■ DR7 

o Las Daraas 



i<r 3 - 



io-" - " 

20 40 60 80 100 120 140 160 180 

r (Mpc/h) 

Figure A4. £o measured from 80 Las Damas mock catalogs and from SDSS 
DR7 data. 




• DR7 

o Las Damas 



10- 3 - 



20 40 60 80 100 120 140 160 180 

r (Mpc/h) 

Figure A5. £2 measured from 80 Las Damas mock catalogs and from SDSS 
DR7 data. 

Bertacca D., Bartolo N., Diaferio A., Matarrese S., JCAP, 0810, 
023, 2008 

Bertacca D., Bartolo N., Matarrese S., Advances in Astronomy, 
904379, 2010 

Bertacca D., Raccanelli A., Piattella O., Pietrobon D., Bartolo N., 
Matarrese S., Giannantonio T., JCAP 1 103, 039, 201 1 

Bertacca D., Maartens R., Raccanelli A., Clarkson C, arXiv: 
1205.5221 

Beynon E., Bacon D. J., Koyama K., MNRAS, 403, 353, 2010 

Bharadwaj S., ApJ, 472, 1, 1996 

Blake C. et al, MNRAS, 415, 2876 B, 201 1 

Brans C. H., arXiv:gr-qc/0506063| 

Cabre A., Gaztanaga E., MNRAS, 393, 1 183, 2009 

Camera S. et al. 2009, MNRAS, 399, 1995 

Camera S. et al. 201 1, MNRAS, 415, 399 

Camera S., Carbone C, Moscardini L., JCAP 1203, 039, 2012 

Carroll S. M., Duvvuri V., Trodden M., Turner M. S., Ph Rv D, 

70d3528C, 2004 
Cole S., Fisher K. B., Weinberg D. H., MNRAS, 275, 515, 1995 
Crocce M., Scoccimarro R., PRD73, 063519, 2006 
Crocce M., Scoccimarro R., PRD77, 023533, 2008 
Dalai N., Dore O., Huterer D., Shirokov A., PhRvD, 77, 123514, 

2008 

de Rham C, Hofmann S., Khoury J., Tolley A. J., Journal of Cos- 
mology and Astro-Particle Physics, 2, 1 1, 2008 



3 




50 100 



r (Mpc\h) 

Figure A6. £0 measured from 80 Las Damas mock catalogs and theoretical 
predictions for ACDM+GR . 



• Mock Measurements 
ACDM 




r (Mpc\h) 

Figure A7. £2 measured from 80 Las Damas mock catalogs and theoretical 
predictions for ACDM+GR . 



Deffayet C, Physics Letters B, 502, 199, 2001 
Desjacques V. and Seljak U., 2010, Class. Quant. Grav., 27, 
124011 

Durrer R., Maartens R., 2008, Gen.Rel.Grav., 40, 301 
Dvali G. R., Gabadadze G., Porrati, M., Phys. Lett., B, 484, 112, 
2000 

Dvali G, Gabadadze G. and Porrati M., Physics Letters B, 485, 
208, 2000 

Dvali G., Hofmann S., Khoury J., Physical Review D, 76, 8, 2007 
Eisenstein D.J., et al., AJ, 122, 2267, 2001 
Eisenstein D. J., Seo H.-J., White M., ApJ, 665, 14, 2007 
Fukugita M., Ichikawa T, Gunn J. E., Doi M., Shimasaku K., 

Schneider D. P., AJ, 111, 1748, 1996 
Gregory R., Kaloper N., Myers R. C, Padilla A., Journal of High 

Energy Physics, 10, 69, 2007 
Gunn J.E., et al, AJ, 1 16, 3040, 1998 
Guzzo L. et al., 2009, Nat., 541, 451 
Hamilton A.J.S., 1992, ApJ, 385, L5 

Hamilton A. J. S., "Linear redshift distortions: A review", in "The 
Evolving Universe", ed. D. Hamilton, pp. 185-275 (Kluwer Aca- 
demic, 1998) |astro-ph/970"8102| 

Hawkins E. et al., MNRAS, 346, 78, 2003 

Hu W., Sawicki I., Phys. Rev. D 76 064004, 2007 

Jackson J.C., MNRAS, 156, 1, 1972 



© 0000 RAS, MNRAS 000, 000-000 



12 Raccanelli et al. 



Jain B„ Zhang P., Phys Rev D, 78, 063503, 2007 
Nesseris S., Perivolaropoulos L., PhRvD..77b3504N, 2008 
Kaiser N., 1987, MNRAS, 227, 1 
Kazin E., et al., 2010, ApJ, 710, 1444 

Koyama K. & Maartens R., Journal of Cosmology and Astro- 
Particle Physics, 1,16, 2006 
Kwan J., Lewis G. R, Under E. V., 2011, |arXiv:1105.1194| 
Landy S.D., Szalay A.S., 1993, ApJ, 412, 64 
Linder E., 2005, PRD, 72, 043529 
Linder E., Astropart. Phys., 29, 336, 2007 

Lombriser L., Hu W., Fang W., Seljak U., arXiv:astro- 

ph/0905.1 112, 2009 
Luty M.A., Porrati M., Rattazzi R. Journal of High Energy 

Physics, 09, 029, 2003 
Marulli E, Bianchi D., Branchini E., Guzzo L., Moscardini L., 

Angulo R. E., 2012, |arXiv:1203 .1002 
Matarrese S., Pietroni, M., JCAP, 0706, 026, 2007 
Matarrese S., Verde L., 2008, ApJ, 677, L77 
Matsubara T., 2004, ApJ, 615, 573 
Matsubara T., PRD, 77, 063530, 2008a 
Matsubara T„ PRD, 78, 083519, 2008b 
Maartens R., Living Rev. Rel., 13, 5, 2010 
McBride C, et al., 201 1, in prep. 

Nicolis A. & Rattazzi R., Journal of High Energy Physics, 6, 59, 
2004 

Okumura T., Matsubara T., Eisenstein D.J., Kayo I., Hikage C, 

Szalay A.S., Schneider D.P, 2008, ApJ, 676, 889 
Papai P., Szapudi I., 2008, MNRAS, 389, 292 
Peacock J. A., Dodds S. J., 1996, MNRAS, 280, L19 
Peacock J.A., et al., 2001, Nature, 410, 169 
Percival W.J., et al., MNRAS, 327, 1297, 2001 
Percival W. J. et al., MNRAS, 385, L78, 2004 
Percival W.J., White M., MNRAS, 393, 297, 2009 
Perlmutter S., et al, 1999, ApJ, 517, 565 
Piattella O., Bertacca D., MPLA, 26, 2277, 201 1 
Pope A.C., et al., ApJ, 607, 655, 2004 

Raccanelli A., Samushia L., Percival W. J., 2010, MNRAS, 409, 
1525 

Reid B. A., et al., MNRAS, 404, 60, 2010 

Reid B.A., White M., 2011, MNRAS, 417, 1913, 
larXiv: 1105 .41651 

Reid B.A., et al., 2012, |arXiv:1203.6641| 
Riess A.G., et al., 1998, AJ, 116, 1009 
Ross A., et al, 2011, |arXiv:1105.2320| 
Ross A., et al, 2012, |arXiv:1203.6499| 

Samushia L., Percival W. J., Raccanelli A., 2012, MNRAS, 420, 
3,2102 

Sanchez A. G., et al., 2012, |arXiv:1203.6616| 
Schmidt E, Physical Review D, 80, 12, 2009 
Scoccimarro R., 2004, PRD, 083007 

Slosar A., Hirata C, Seljak U., Ho S., Padmanabhan N., 2008, 
JCAP 08, 031 

Song Y.-S., Sawicki I. and Hu W., Physical Review D, 75, 6, 2007 

Song Y.-S., Koyama K., JCAP, 01, 048, S 2009 

Song Y.-S., Percival W. J., JCAP, 0910, 004, 2009 

Szalay A.S., Matsubara T„ Landy S.D., 1998, ApJ, 498, LI 

Szapudi I., 2004, PhRvD, 70, 083536 

Taruya A., Nishimichi T., Saito S., 2010, PRD, 063522 

Tegmark M. et al., PRD, 74, 123507, 2006 

Thomas S. A., Abdalla F. B., Lahav O., 2010, |arXiv:1012.2272| 

Tsujikawa S., Lect. Notes Phys., 800, 99-145, 2010 

Varshalovich D. A., Moskalev A. N., & Khershonski V. K, 1988, 



Quantum Theory of Angular Momentum (Singapore: World Sci- 
entific) 

Wang Y, 2012, |arXiv: 1202.5238| 

Xia J.-Q., Bonaldi A., Baccigalupi C, De Zotti G., Matarrese S., 

Verde L., Viel M., 2010, JCAP, 1008, 013 
Yamamoto K., Sato T, Hiitsi G., Progress of Theoretical Physics, 

120, 3, 609, 2008 
Yamamoto K, Nakamura G., Hiitsi G., Narikawa T, Sato T, 

Physical Review D, 81, 10, 103517, 2010 
Zaroubi S., Hoffman Y, 1993, |arXjv:astro-ph/93110l3l 
Zehavi I. et al., ApJ, 621, 22, 2005 



© 0000 RAS, MNRAS 000, 000-000 



